One-loop calculation of mass dependent O(a) improvement coefficients 
for the relativistic heavy quarks on the lattice 
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We carry out the one-loop calculation of mass dependent O(a) improvement coefficients in the relativistic heavy 
quark action recently proposed, employing the ordinary perturbation theory with the fictitious gluon mass as an 
infrared regulator. We also determine renormalization factors and improvement coefficients for the axial- vector 
current at the one-loop level. It is shown that the improvement coefficients are infrared finite at the one-loop level 
if and only if the improvement coefficients in the action are properly tuned at the tree level. 



1. INTRODUCTION 

While, in principle, lattice QCD allows a pre- 
cise determination of weak matrix elements asso- 
ciated with the B and D mesons, c- and 6-quarks 
are too heavy to treat directly in current numer- 
ical simulations due to large 0(niQa) errors. 

Recently, a new relativistic approach to control 
niQa errors was proposed from the view point of 
the on-shell 0(a) improvement program assuming 
rriQ ;§> Ancn [T12] . The action is given by 

S q = a 4 Y<xQ( x ) [loD + v r )iD l + mo 

~ar t /2Dl - ar s /2D? - aigc B / Aa^F^ 

-aigc E /2a ai Foi]q{x), (1) 

where we are allowed to choose r t = 1 and other 
four parameters is, r s , ce and cb are analytic 
functions of mQa and the gauge coupling g. Once 
a nonperturbative adjustment of four parameters 
v, r s , ce and cb is achieved, the remaining cutoff 
effects are reduced to be ©((oAqcd) 2 )- In this re- 
port we present a perturbativc determination of 
ce and cb up to the one- loop level (see Ref. for 
details). In addition we calculate the renormal- 
ization factors and the 0{a) improvement coeffi- 
cients for the axial-vector current at the one-loop 
level. 



2. DETERMINATION OF c E AND c B AT 
THE ONE-LOOP LEVEL 

We determine c E — effl + g 2 c^ and cb — 
c^P + g 2 c B 1 ^ from the on-shell quark-quark scat- 
tering amplitude. In the massless case, csw is 
successfully determined up to the one-loop level 
employing the conventional perturbation theory 
with the fictitious gluon mass as an infrared reg- 
ulator [3| . We extend this calculation to the mas- 
sive case. Since the tree level determination of eg 
and cs in the massive case is already performed 
in Ref. pQ, next step is to calculate six types of 
one-loop diagrams depicted in Fig. Q 

We first consider to calculate . Without 
the space-time symmetry the general form of the 
off-shell vertex function at the one-loop level is 
written as 

Afc 1 ' (p, 1, m ) 

= IkF* + 7fc {^ 2 fc + tfoF^} + {4F% + 4 F*} lk 

+4iki>F* + 4 lk hF$ + 4o7krfFg 

+ ( Pk + q k ) [H k ; + i>H k + 4H* + 4tfH*] (2) 
+ (Pk - qk) [G\ + $G k 2 + 4G\ + 4iG k ^ + 0(a 2 ), 
where A fc = A^ 0) + g 2 A^ ] + 0(g 4 ) and j = 

J2 3 a =oPaJa, 4 = X)c,=0?a7a> $0 = Po70, 4o = 

9o7o- P and q are incoming and outgoing quark 




Figure 1. One-loop quark-gluon vertex diagrams. 



momenta, respectively. The coefficients F*, G\ 
and H\ are functions of p 2 , q 2 , p - q and m. Sand- 
wiching Al 1 (pi q, m) by the on-shell quark states 
u(p) and u(q), which satisfy $u(p) — im p u(p) and 
u(q)4 — im p u(q), the matrix element is reduced 
to 

u(q)A^ ) (p, q,m)u(p) 

= u{q) lk u{p) {f* + im p 2F$ - m 2 p Fl} 

+ {Pk + q k )u{q)u{p) {iff + im p 2H^ 

+ 0(a 2 ), (3) 

where we use the relations among F* , G\ and Hf 
constrained by the charge conjugation symmetry. 

The relevant term for the determination of eg 
is the last one in eq.Q, whose coefficient can be 
extracted by setting p = p* = (po = im p ,pi = 0) 
and q = q* = (qo = im pi qi = 0) in cq.JSJ: 

[Hi + im p (H$ + H$) - m 2 p H k 4 ][Z , a=a . 
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1) (4) 



We should remark that the second term 
in eq.@ contains both the lattice artifact of 
0(pka, qka) and the physical contribution of 
0(jJk/m, qk/m). The parameter cb is determined 
to eliminate the lattice artifact of 0(p k a, qka): 

C * = [H} + imJH* - m 2 H*] la " . . 
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Fi gure 2. and c^ 1 with the plaquette gauge 
action as a function of pole mass m p °K 



-Zf) [iff + im p 2H* - m 2 p Hf\ 



cont 

p=p* ,q=q* 



The infrared behavior of eq.Q} is investigated 
by expanding the inner loop momentum in Fig. ^ 
around zero. The infrared divergence is found to 
be 



2N r + 2 l [ ° B 



r(P) 
•q 



where L = -1/(16tt 2 ) In |A 2 a 2 | with A the ficti- 
tious gluon mass. If the tree level improvement 
coefficients are properly tuned as — ri°^ , we 
are left with -{N c /2){zf ] /m p 0) )L, which is ex- 
actly the same as the infrared divergence in the 
continuum theory with the correct normalization 



factor Zq°\ 

In the same way, 
the infrared divergence. Results for Ca J and c 



e is determined free from 
(i) „„j Ji) 



with the plaquette action are plotted in Fig.[5| As 
expected c^ 1 ' and Cjp approaches to the massless 
value for cswEO as m^a vanishes. 



0(a) IMPROVEMENT OF 
AXIAL- VECTOR CURRENT 



THE 



We consider the heavy-light axial-vector cur- 
rent Ap. Without space-time symmetry, the 
renormalized operator A^ tt,R {x) with the 0(a) 
improvement is written as 

A^> R {x)=z¥*[q{x) W Q{x) 
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Figure 3. in heavy-light system with the pla- 
quette gauge action as a function of heavy quark 
pole mass /rig . 



-9 2c X d t il( x h5Q(x)} - 9 2 c A ^ {q{x)j 5 Q(x)} 
'9 2c AAdzq(x)}^f,l5Q(x) 



2 cf u q(x)'Y tl 'Ys'Yi{diQ(x)} 



(5) 



where Z A U and 



depend on the quark 
masses ttiq and m q . and d~ are defined as 



9, 



5 P and = <9 M - 



d„. With the aid of 



equation of motion we can choose c 



h — L 



0. 



To extract Z 1 ^ and 



c t' < a ^i qx the one-loop 
level from the general form of the vertex func- 
tion, we have to use both decay and scattering 
processes. Except for the physical contribution, 



-,H,L) 



all parameters are infrared finite, once , 



„(°) 



...(") 



and are properly tuned. In Figs. 0] 



<+,-,H,L) 



we show numerical results for and 
in heavy-light system. The former is defined by 

where we take 



^latt /^MS _ 



4 0) (i-5 2 A, 



J A I "A 

m q = for the light quark 



We observe antici- 
pated features that c^' H ' L ' vanishes as the heavy 

quark mass decreases, while A^ and c\ become 
close to their corresponding massless values ob- 
tained in Refs. |4I5| . 
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Figure 4. c^' , c^' ' H ' L ^ in heavy-light system 
with the plaquette gauge action 
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